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The  flow  field  generated  by  a  zero-net  mass-flux  (ZNMF)  actuator  is  investigated  via 
both  numerical  simulations  and  experiments  to  augment  the  current  understanding  of  the 
flow  physics  of  the  orifice.  The  results  aid  in  improving  the  accuracy  of  low  dimensional 
lumped  element  ZNMF  models  suitable  for  design.  Dimensional  analysis  yields  a  number  of 
key  parameters  that  govern  the  characteristics  of  this  flow.  Among  them  for  a  sharp 
rectangular  slot  or  circular  orifice  are  the  Reynolds  number,  the  dimensionless  stroke 
length,  and  the  orifice  height-to-diameter  ratio.  Variation  of  these  parameters  shows  that 
the  flow  field  differs  appreciably  from  the  exact  linear  solution  of  pipe  flow  driven  by  an 
oscillatory  pressure  gradient.  In  particular,  depending  on  the  stroke  length  and  the  orifice 
geometry,  the  pressure  drop  in  the  orifice  may  be  dominated  by  nonlinear  “minor”  losses 
due  to  entrance/exit  effects,  or  linear  “major”  losses  associated  with  the  presence  of  a 
nominally  fully-developed  region  in  the  central  region  of  the  orifiee/slot. 


I.  Introduction 

ZERO-net  mass-flux  (ZNMF)  actuators,  also  known  as  synthetic  jets,  have  emerged  as  versatile  actuators  in 
various  applications  (e.g.,  thrust  vectoring  of  jets,1  active  control  of  separation2'34,5  and  turbulence  in  boundary 
layers  ’  ).  ZNMF  actuators  exhausting  into  a  quiescent  medium  have  been  studied  both  experimentally  ’  ’  and 
numerically,11,12  while  other  studies  have  focused  on  the  interaction  with  an  external  boundary  layer.13,14,15,16 
However,  numerous  unresolved  issues  remain  concerning  the  fundamental  governing  physics  of  these  devices, 
effectively  hindering  their  modeling,  design  and  optimization.  For  example,  the  unsteady  flow  in  the  orifice  or  slot 
plays  a  large  role  in  determining  the  actuator  performance.  While  numerous  parametric  studies  have  examined 
various  orifice  geometry  and  flow  conditions,  a  clear  understanding  of  the  loss  mechanisms  is  still  lacking.  Detailed 
numerical  simulations  (and  companion  experiments)  can  be  used  to  elucidate  the  underlying  physics  but  are  not 
practical  as  a  design  tool.  Instead,  accurate  low  dimensional  models  are  required  to  facilitate  the  effective  design  of 
ZNMF  actuators  for  specific  applications. 

Fig.l  shows  a  simple  schematic  of  a  typical  ZNMF  actuator.  Fluctuating  flow  in  the  slot  is  induced  by  the 
oscillation  of  the  driver  mounted  in  a  cavity  of  volume  V  using,  for  example,  a  piezoelectric  diaphragm  or  an 
electromagnetic  piston,  etc.  The  diaphragm  deflection  is  characterized  by  its  volume  displacement  AV  as  a 
function  of  the  angular  frequency  co .  The  orifice  geometry  is  described  by  its  height  h  and  diameter  d  for  a 
sharp-edged  axisymmetric  orifice  and  similarly  for  a  two-dimensional  slot. 
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With  regards  to  the  nature  of  the  unsteady  flow  in  the  orifice 
and  the  associated  losses,  it  is  instructive  to  first  consider  the 
simpler  case  of  steady  flow  through  a  pipe.  In  this  case,  losses 
arise  due  to  different  mechanisms.  In  any  undergraduate  fluid 
mechanics  textbook,  these  losses  are  characterized  as  “major” 
losses  in  the  fully  developed  flow  region  and  “minor”  losses 
associated  with  entrance  and  exit  effects,  etc.  For  laminar  flow, 
the  pressure  drop  A p  in  the  fully-developed  region  is  linearly 
proportional  to  the  volume  flow  rate  Q  or  average  spatial 
velocity  uavg  ,  while  the  nonlinear  minor  pressure  losses  are 
proportional  to  the  dynamic  pressure  0.5 pu2  .  Similarly,  for 

the  case  of  unsteady ,  laminar,  fully-developed,  flow  driven  by 
an  oscillatory  pressure  gradient,  the  complex  flow  impedance, 

Ap'/Q',  can  be  determined  analytically  and  decomposed  into  linear  resistance  and  reactive  components.17 
Unfortunately,  no  such  solution  is  available  for  the  nonlinear,  and  perhaps  dominant,  losses  associated  with  entrance 
and  exit  effects.  The  lumped  element  model  of  McCormick18  uses  the  limiting  steady  flow  solution  for  Ap'/Q' , 
while  that  of  Gallas  et  al.19  uses  the  frequency-dependent  analytical  solution  for  the  linear  resistance.  However,  both 
models  suffer  from  an  overly  simplistic  empirical  equation  for  the  nonlinear  resistance  component. 

An  alternative  approach  to  characterize  or  “calibrate”  an  oscillatory  fluidic  actuator  uses  the  empirical 
observation  that  the  cavity  pressure  fluctuation  A p'  is  linearly  proportional  to  the  centerline  exit  velocity  fluctuation 

u  at  low  forcing  levels  and  oc  u'2  (i.e.,  nonlinear)  at  sufficiently  high  forcing  levels.20'21  However,  the  relative 
importance  and  scaling  of  the  linear  and  nonlinear  components  versus  the  governing  dimensionless  parameters  is 
unknown  and  remains  a  critical  obstacle  for  designers  of  ZNMF  actuators. 

This  paper  presents  the  progress  in  an  ongoing  collaborative  effort  focused  on  the  development  of  low¬ 
dimensional  models  for  ZNMF  actuators.  While  the  ultimate  goal  is  to  develop  models  suitable  for  use  in  boundary 
layer  flow  control  in  which  the  ZNMF  actuator  interacts  with  a  grazing  flow,  the  current  paper  focuses  on 
experimental  and  computational  efforts  to  first  understand  and  model  the  oscillatory  orifice  flow  in  the  simpler  case 
without  a  grazing  flow.  As  such,  similar  to  the  recent  work  by  Smith  &  Swift  who  experimentally  studied  the  losses 
in  an  oscillatory  flow  through  a  rounded  slot,2’  this  paper  focuses  on  the  more  commonly  employed  sharp-edged 
orifice  or  slot  of  a  ZNMF  actuator.  In  particular,  experimental  data  and  numerical  simulations  are  used  to  explore 
the  flow  behavior  and  develop  components  for  a  low-order  lumped  element  model  of  a  ZNMF  actuator.  The  effects 
of  various  governing  dimensionless  parameters  are  examined.  The  paper  is  organized  as  follows.  Following  a 
dimensional  analysis,  the  details  of  the  experimental  and  computational  methods  are  presented.  The  remainder  of 
the  paper  focuses  on  analysis  of  the  resulting  data  to  elucidate  the  flow  physics  and  to  develop  a  model  for  nonlinear 
resistance.  Finally,  preliminary  conclusions  are  offered,  and  future  work  is  discussed. 

II.  Nondimensional  Parameters 

The  geometry  of  the  flow  under  consideration  is  shown  in  Fig.  1.  The  governing  dimensionless  parameters  for  a 
ZNMF  actuator  include  a  Reynolds  number  and  either  a  normalized  stroke  length  or  a  Strouhal  number.  Various 
velocity  scales  are  employed  in  the  literature,  but  the  present  work  employs  the  spatial-averaged,  time-averaged  exit 
velocity  during  the  expulsion  portion  of  the  cycle: 


x 


Figure  1.  Schematic  of  a  typical  ZNMF 
actuator. 


u)  =~L\Tu^y)dtdA 
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where  T  =  1/ /  =  2njw  is  the  period  of  the  cycle,  An  is  the  orifice  area,  and  u  (t)  is  defined  as  the  spatial- 
averaged,  axial  component  of  the  velocity  at  the  exit  plane.  The  stroke  length  L0 ,  which  is  a  measure  of  the  distance 
that  a  fluid  particle  travels  during  the  expulsion  part  of  the  cycle,  is  thus  defined  as 

Lo=\Tu^)dt  =  tJJT2-  (2) 
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Note  the  equivalence  between  the  different  dimensionless  parameters.  The  Strouhal  number  is  defined  as 
St  =  codj  Uj  ;  thus  the  following  relationship  exists  between  the  Strouhal  number,  Stokes  number  S  =  -Jcod^Jv  , 
Reynolds  number  Re7  =  Ujd  jv  ,  and  the  dimensionless  stroke  length  L0/d  : 


Re./  _  Ljd 
S2 


A  dimensional  analysis  of  a  ZNMF  actuator  with  a  sharp-edged  rectangular  2-D  slot  or  orifice  using  a  generic 
driver  and  exhausting  in  a  quiescent  medium  reveals  the  following  functional  form: 


(»,Re,)  =  / 


co  h  co 

'T  ’  ’ 

coH  d  cod 


(4) 


The  two  parameters  on  the  LHS  are  both  primarily  dependent  on  Uj .  The  terms  on  the  RHS  are  explained  as 
follows:  The  ratio  of  the  driving  frequency  to  the  Helmholtz  frequency  is  a  measure  of  the  compressibility 

of  the  flow  inside  the  cavity.  The  flow  in  the  cavity  is  essentially  incompressible  provided  that  co/coH  <SC  1 .  Hence, 
incompressible  CFD  models  are  only  relevant  for  this  case.  Note  that  coH  =  c0JAjh'\/  is  based  on  the  isentropic 
sound  speed  c0 ,  the  neck  area  An ,  the  effective  height  of  the  neck  h' ,  and  cavity  volume  V  .  Hence,  the  condition 
co/wH  <£.  1  is  satisfied  for  a  sufficiently  small  cavity  volume.  Next,  co/cod  is  the  ratio  of  the  drive  frequency  to  the 
natural  frequency  of  the  driver,  and  A \//d2  is  the  dimensionless  displacement  of  the  driver.  The  parameter 
kd  =  IndjA.  is  related  to  the  ratio  of  the  orifice  diameter  to  the  acoustic  wavelength.  The  Stokes  number  S  can  be 
interpreted  as  the  ratio  of  the  orifice  diameter  to  the  unsteady  boundary  layer  thickness  in  the  orifice  Jv/co  .  Eq.  (4) 
forms  the  basis  for  the  selection  of  our  test  matrix,  which  for  the  purposes  of  studying  the  oscillatory  flow  in  the 
orifice,  consists  of  h/d  ,  Rey ,  and  S  (or,  by  virtue  of  Eq.  (3),  L0/d  or  St ). 

III.  Experimental  Setup 

Fig.  2  shows  a  schematic  of  the  experimental  setup. 

A  piezoelectric  diaphragm  is  actuated  using  an  Agilent 
33120A  function  generator  with  an  amplifier.  Cavity 
pressure  measurements  are  obtained  using  a  flush- 
mounted  Briiel  and  Kjasr  (B&K)  1/8"  diameter 

condenser  type  microphone  (Model  4138)  powered  by  a 
B&K  2670  pre-amplifier  and  a  B&K  2804  power 
supply.  In  addition,  the  centerline  deflection  of  the 
diaphragm  is  measured  using  a  laser  displacement 
sensor  (Micro-Epsilon  Model  ILD2000-10).  Both 
signals  are  acquired  phase-locked  to  the  drive  signal 
and  sampled  simultaneously  with  a  dynamic  signal 
analyzer  (National  Instillments  Model  4552). 

Velocity  measurements  of  the  tlowfield  emanating 
from  the  ZNMF  orifice  are  acquired  using  Laser 
Doppler  Anemometry  (LDA),  the  details  of  which  are 
listed  in  Table  1.  The  synthetic  jet  actuator  is  mounted 
to  a  three-axis  traverse  with  sub-micron  spatial 
resolution.  The  488  and  514.5  nm  lines  of  an  argon-ion 
laser  are  used  to  obtain  coincident,  two-component 
velocity  measurements  using  a  Dantec  FiberFlow 
system.  As  shown  in  Fig.  3,  a  novel  three-beam  optical  Figure  2.  Schematic  of  the  experimental  setup  for  off- 
combiner  configuration  is  used  to  facilitate  velocity  axis,  two-component  LDA  measurements. 
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Table  2.  LDA  measurement  details. 


Probe 


Figure  3.  LDA  3-beam  optical  configuration. 


Property: 

LDA  1 

LDA  2 

Wavelength  (nm) 

514.5 

488 

Focal  length  (mm) 

120 

120 

Beam  diameter  (mm) 

1.35 

1.35 

Beam  spacing  (mm) 

26.87 

26.87 

Number  of  fringes 

25 

25 

Fringe  spacing  (pm) 

2.3121 

2.1929 

Beam  half-angle  (deg) 

6.39 

6.39 

Probe  volume  -  dx  (mm) 

0.0586 

0.0556 

Probe  volume  -  dy  (mm) 

0.0582 

0.0552 

Probe  volume  -  dz  (mm) 

0.5233 

0.4963 

Table  1.  Synthetic  jet  actuator  details. 


Property: 

Value: 

Cavity  volume  V  (m3) 

7.11xl0"6 

Orifice  diameter  d  (mm) 

2.00 

Orifice  thickness  h  (mm) 

1.65 

measurements  at  the  exit  plane  of  the  synthetic  jet 
actuator  (see  Table  2).  A  direction  cosine 
transformation  is  then  applied  to  the  acquired 
velocity  components  LDA  1  and  LDA  2  to  extract 
the  axial  and  radial  velocity  components. 

A  200  mm  micro  lens  and  bellows  extender 
collects  lights  at  90°  off-axis  in  order  to  improve  the 
spatial  resolution.  Scattered  light  from  the  probe 
volume  is  focused  and  passed  through  a  100  /jm  diameter  pinhole  aperture.  The  resulting  field  of  view  was  imaged 
using  a  micro-ruler  and  found  to  be  approximately  10  //in.  indicating  that  the  effective  length  of  the  probe  volume 
dz  has  been  reduced  by  over  an  order  of  magnitude  from  that  listed  in  Table  1.  After  the  pinhole,  a  color  separator 
splits  the  514.5  nm  and  488  nm  wavelengths  and  transmits  the  light  to  two  separate  photomultiplier  tubes  (PMTs) 
for  subsequent  FFT  analysis. 

Phase-locked  velocity  measurements  are  acquired  at 
the  centerline  of  the  jet  exit  with  a  resolution  of  10°  and 
approximately  200  velocity  values  are  acquired  at  each 
phase  angle.  In  addition,  the  synthetic  jet  device  is 
traversed  to  acquire  phase-locked  velocity  profdes 
across  the  surface  of  the  orifice  with  a  spatial  resolution 
of  0. 1  mm  for  select  cases. 

Once  the  velocity  profiles  at  the  surface  of  the 
orifice  are  acquired,  the  phase-locked  profiles  are 
spatially  integrated  to  determine  the  periodic  volume 
flow  rate.  The  velocity  profiles  are  phase-locked  to  the 
other  measured  signals  via  the  TTL  trigger  signal 
available  on  the  function  generator.  In  the  data 
presented  in  this  paper,  zero  phase  is  defined  as  the 
point  at  which  the  volume  flow  rate  is  zero  with 
positive  slope. 

Fig.  4  illustrates  a  set  of  typical  phase-locked  axial 
velocity  profiles  during  four  different  phases  separated 
by  90°  in  the  cycle,  corresponding  to  maximum 

expulsion,  maximum  ingestion,  and  the  two  phases  half  way  between.  The  error  bars  represent  an  estimate  of  the 
95%  confidence  interval  for  each  velocity  measurement  and  are  obtained  using  a  perturbation  technique23  that  yields 
the  same  nominal  values  of  uncertainty  as  a  standard  Monte  Carlo  technique  but  with  significantly  less 
computational  time.  This  method  is  employed  to  estimate  the  uncertainty  in  the  volume  flow  rate.  The  circular 
orifice  was  observed  under  a  microscope  and  found  to  have  a  diameter  d  =  2  ±  0.06  mm.  The  95%  confidence 
intervals  estimates  of  Qt  and  d  ,  in  turn,  are  used  to  estimate  the  uncertainty  in  the  Reynolds  number,  which  is  in 
the  range  of  2-10%. 
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Figure  4.  Phase-locked  axial  velocity  profiles  at 
x/d=0.07,  S=8,  Re,  =45.5. 
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The  locus  of  the  positive  values  of  the  volume  flow  rate  are  integrated  to  give  the  average  volume  flow  rate 
during  the  expulsion  part  of  the  cycle,  Qj  ,  which  is  related  to  the  average  velocity  by  Qf  =  U.nd2  j  A  .  Thus,  the 
Reynolds  number  is  given  by 


40: 

Re7  =-=i- 

1  7T\’d 


IV.  Numerical  Methodology 

An  in-house,  cell-centered  based  Cartesian  grid  solver  is  used  for  simulating  the  flow  generated  by  a  ZNMF 
actuator.  The  incompressible  Navier-Stokes  equations  are  used  since  <K  1  and  are  non-dimensionalized  with 

appropriate  length  and  velocity  scales  and  discretized  using  a  cell-centered,  collocated  (non-staggered)  arrangement 
of  the  primitive  variables  velocities  u  and  pressure  p  .  In  addition  to  the  cell-center  velocities  ( u  ),  the  face-center 


velocities  U  are  computed.  Similar  to  a  fully  staggered  arrangement,  only  the  component  normal  to  the  cell-face  is 
calculated  and  stored.  The  face -center  velocity  is  used  for  computing  the  volume  flux  from  each  cell.  The 
advantage  of  separately  computing  the  face -center  velocities  was  initially  proposed  by  Zang  et  al.~4  and  discussed  in 
the  context  of  the  current  method  in  Ye  et  al.25  A  second- 


order  Adams-Bashforth  scheme  is  employed  for  the 
convective  terms,  while  the  diffusion  terms  are  discretized 
using  an  implicit  Crank-Nicolson  scheme  which  eliminates 
the  viscous  stability  constraint.  The  pressure  Poisson 
equation  is  solved  with  a  Krylov-based  approach. 

A  multi-dimensional  ghost-cell  methodology  is  used  to 
incorporate  the  effect  of  the  immersed  boundary  on  the 
flow.26  The  schematic  in  Fig.  5  (a)  shows  a  solid  body  with 
a  curved  boundary  moving  through  a  fluid  and  the 
categorization  of  the  nodes  for  the  purpose  of  implementing 
the  ghost  cell  methodology.  For  moving  boundaries,  the 
general  framework  can  be  considered  as  Eulerian- 
Lagrangian,  wherein  the  immersed  boundaries  are  explicitly 
tracked  as  surfaces  in  a  Lagrangian  mode,  while  the  flow 
computations  are  performed  on  a  fixed  Eulerian  mesh. 
Hence,  we  identify  cells  that  are  just  inside  the  immersed 
boundaries  as  “ghost  cells.”  The  discrete  equations  for  these 
cells  are  then  formulated  as  to  satisfy  the  imposed  boundary 
condition  on  the  nearby  flow  boundary  to  second-order 
accuracy.  These  equations  are  then  solved  in  a  fully  coupled 
manner  with  the  governing  flow  equations  of  the  regular 
fluid  cells.  Care  has  been  taken  to  ensure  that  the  equations 
for  the  ghost  cells  satisfy  local  and  global  mass  conservation 
constraints  as  well  as  pressure-velocity  compatibility 
relations.  The  solver  has  been  designed  to  take  geometrical 
input  from  conventional  CAD  programs.  The  code  has  been 
well  validated  by  comparisons  against  established 
experimental  and  computational  data.27  The  solver  is 
currently  being  used  to  carry  out  parametric  studies  of 
synthetic  jets  in  both  quiescent  and  grazing  flow  fields.14 
The  typical  2-D  setup  for  a  rectangular  synthetic  jet  in 
quiescent  flow  field  is  as  shown  in  Fig.  5  (b).  For  the 
numerical  simulations  described  here,  a  pulsatile  boundary 
condition  is  provided  at  the  bottom  of  the  cavity. 


Figure  5.  (a)  Schematic  of  the  sharp-interface 
method  on  a  fixed  Cartesian  mesh  and  (b) 
synthetic  jet  computational  boundary  conditions. 
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V.  Lumped  Element  Modeling 

The  low  dimensional  modeling  approach  adopted  is  based  on  lumped  element  modeling  (LEM).18 19  In  LEM,  the 
individual  components  of  a  synthetic  jet  are  modeled  as  elements  of  an  equivalent  electrical  circuit  using  conjugate 
power  variables  (i.e.,  power  =  generalized  “flow”  x  generalized  “effort”  variables).  Fig.  6  shows  an  equivalent 
circuit  representation  of  a  piezoelectric-driven  synthetic  jet  actuator,  where  the  lumped  parameters  represent 
generalized  energy  storage  elements  (i.e.,  capacitors  and  inductors)  and  dissipative  elements  (i.e.,  resistors).  Model 
parameter  estimation  techniques,  assumptions,  and  limitations  are  discussed  in  Gallas  et  al.|Q  The  frequency 


Figure  6.  Equivalent  circuit  model  of  a  piezoelectric- 
driven  synthetic  jet  actuator. 


response  function  of  the  circuit  is  derived  to  obtain  an 
expression  for  Q}  jVac  ,  the  jet  volume  flow  rate  during 

the  expulsion  part  of  the  cycle  per  applied  voltage. 
LEM  provides  a  compact  nonlinear  analytical  model 
and  valuable  physical  insight  into  the  dependence  of 
the  device  behavior  on  geometry  and  material 
properties.  However,  the  primary  difficulty  with  this 
method  is  the  lack  of  sufficient  physics  in  the  orifice 
impedance  models,  embodied  as  linear  and  nonlinear 
resistors  and  masses,  as  shown  in  Fig.  6. 


Figure  7.  Volume  velocity  divider  circuit,  valid 
for  a  generic  ZNMF  actuator. 


With  slight  modifications,  the  model  is  applicable  to  any  type  of  ZNMF  actuator  (electromagnetic,  etc.)  and  for 
two-dimensional  slots  and  axisymmetric  orifices.  Fig.  7  shows  a  current  divider  circuit  that  is  applicable  to  a 
generic  driver  provided  that  the  motion  (i.e.,  the  volume  velocity  Qd  )  of  the  driver  is  known. 

The  sinusoidal  volume  velocity  of  the  driver  is  Qd  =  /V/jAV sin  (on ) ,  where  AV  is  the  net  volume  displaced  by 
the  driver  motion.  Qd  is  split  into  two  paths,  and  the  flow  takes  the  path  of  least  impedance.  Essentially,  the  driver 
can  compress  or  expand  the  fluid  in  the  cavity  ((Jd  -  Q.  )  and  expel  or  ingest  fluid  through  the  orifice  ( Qt ) .  The 
governing  equation  is  given  by 


Qj  = 


Q,lZa. 


7  +7 

^  aO  T  ^ aC 


(6) 


where  ZaC  =  1/  jcoCaC  is  the  acoustic  impedance  of  the  cavity,  and  the  acoustic  impedance  of  the  orifice  is  given  by 
ZaQ  ~  Ron  +  R-aOni  +  J 00  iM aN  +  M aRad )  •  ^  desired,  elementary  circuit  laws  can  be  used  to  derive  a  nonlinear  2nd- 
order  ODE  that  describes  the  flow  in  the  orifice. 

The  nonlinear  resistance  Ra0nl  is  associated  with  entrance  and  exit  losses  and  is  usually  modeled  as  a  generalized 
Bernoulli  flow  meter  like  a  minor  “dump”  loss18 

R,m=*pIQj  =  q-5KdpQjIaI>  (?) 
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where  the  dimensionless  loss  coefficient  KD  =  Ap j (0 .5 pu2avg)  was  assumed  to  be  unity  in  ref.  18  and  19.  In  the 
more  general  case,  KD  is  a  function  of  geometry,  Reynolds  number,  and  in  the  case  of  oscillatory  flow,  Stokes 
number  (or  dimensionless  stroke  length  or  Strouhal  number).  Furthermore,  if  the  geometry  is  not  symmetric,  then 
the  loss  coefficient  may  be  different  for  the  inlet  and  exit  portions  of  the  cycle.  Hence,  a  detailed  analysis  on  the 
loss  coefficient  for  various  orifice  shapes,  such  as  that  of  Smith  and  Swift22  for  a  rounded  edge,  should  yield  a  more 
accurate  expression  in  terms  of  modeling  the  associated  nonlinear  resistance.  The  results  presented  below 
summarize  efforts  at  assessing  (1)  the  accuracy  of  Stokes  oscillating  pressure  gradient  solution  to  estimate  the  linear 
impedance  and  (2)  the  variation  of  KD  as  a  function  of  Reynolds  number  and  Stokes  number. 

VI.  Results 

Table  3  shows  the  test  matrix  employed  to  conduct  the  experimental  and  numerical  study  of  the  oscillatory  flow 
through  a  ZNMF  actuator  orifice  by  varying  h/d  ,  Re;  ,  and  S  .  (In  addition,  also  listed  for  convenience  are  /?/ L0 

and  a>la>H  ,  as  well  as  whether  a  jet  is 
formed  or  not  from  the  criterion  based 
on  the  ratio  R Zj / S2  28’  29).  Cases  1-4 

are  experiments,  while  Cases  5-7  are 
computations.  It  should  be  noted  that  10 
additional  experimental  cases  similar  to 
cases  1  and  2  have  been  studied  with 
h/d  =  0.82  and  5=8  by  varying  Re;. 

from  9  to  94.  Similarly,  16  additional 
experimental  cases  similar  to  cases  3  and 

1 

0.8 
0.6 
0.4 

§  0.2 

a 

IS  o 

N 

|  -0.2 

o 

^  -0.4 
-0.6 
-0.8 
-1 


Phase  (deg)  Phase  (deg) 


flow  rate,  (o)  cavity  pressure,  (U)  centerline  driver  velocity. 
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Table  3.  Test  Matrix. 


Case 

Jet  formed 

h/d 

Re, 

5 

h/L o 

co/(oH 

1  (exp.) 

yes 

0.82 

45 

8 

0.37 

0.05 

2  (exp.) 

yes 

0.82 

94 

8 

0.18 

0.05 

3  (exp.) 

no 

0.82 

185 

36 

1.82 

0.93 

4  (exp.) 

yes 

0.82 

1136 

36 

0.29 

0.93 

5  (num.) 

no 

1 

262 

25 

0.76 

0.13 

6  (num.) 

no 

2 

262 

25 

1.52 

0.09 

7  (num.) 

yes 

0.68 

262 

10 

0.08 

0.18 

4  have  been  studied  with  h/d  =0.82  and  S  =36  by  varying  Re;  from  185  to  1875.  However,  only  the  phase  locked 

centerline  velocity  was  acquired  for  these  cases,  while  for  the  first  four  cases  listed  in  Table  3,  velocity  profiles  were 
acquired  across  the  orifice  exit. 

Volume  flow  rate,  cavity  pressure,  and  centerline  driver  velocity  are  shown  phase  locked  in  Fig.  8  for  Cases  1 
through  4,  where  zero  phase  corresponds  to  zero  volume  flow  rate  with  a  positive  slope.  The  measured  driver 
displacement  has  been  differentiated  to  obtain  the  driver  velocity.  The  interesting  point  to  note  is  that,  for  Cases  1 
and  2  with  S  =8  and  a>lwH  =0.05,  the  volume  flow  rate  and  centerline  driver  velocity  are  nearly  in  phase,  indicating 


that  the  flow  is  incompressible.  In  contrast,  for  Cases  3  and  4 
rate  is  not  in  phase  with  the  driver  velocity,  ostensibly 
due  to  compressibility  effects  in  the  cavity. 

The  cavity  pressure  fluctuation  deviates  from  a 
sinusoid  for  Cases  1  and  2,  indicating  significant 
nonlinearity  despite  the  low  Re  .  In  contrast,  in  Cases  3 

and  4,  the  pressure  oscillations  are  nearly  perfect 
sinusoids,  indicating  linear  oscillations  despite  the 
significantly  larger  Re; .  A  possible  explanation  for  the 

observed  behavior  is  discussed  further  below. 

Although  not  shown  here,  a  spectrum  analysis  has 
been  performed  on  the  pressure  traces  in  order  to 
characterize  the  dominant  features  of  the  distorted  time 
signals.  It  was  found  that  only  subharmonics  were 
present,  with  the  total  power  being  dominated  by  the  one 
due  to  the  fundamental  frequency.  Furthermore,  in  order 
to  ascertain  the  resistance  and  reactance  terms  in  the 
orifice,  both  the  amplitude  and  relative  phase  between 
the  pressure  fluctuations  and  volume  flow  rate  at  the 
drive  frequency  are  required.  Fig.  8  also  plots  the 
waveforms  for  both  volume  flow  rate  and  pressure 
determined  from  the  first  three  terms  of  a  Fourier 


S  =36  and  =0.93,  the  orifice  volume  flow 
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Figure  9.  Pressure  fluctuation  normalized  by  dynamic 
pressure  based  on  centerline  velocity  vs.  L0/h:  (*)  S=8, 
9< Re;  <90;  (o)  S=36, 190<Rey  <1875;  (▲)  Case  5;  (■) 
Case  6;  (♦)  Case  7. 


decomposition  via  least  squares  estimation.  Despite  the 
significant  deviation  from  a  sinusoidal 
profile  in  the  pressure  fluctuation  for 
Case  2  for  example,  the  Fourier 
decomposition  provides  a  good 
representation.  The  relative  phase 
between  the  fundamentals  of  the 
pressure  and  volume  flow  rate  can  then 
be  extracted  at  the  fundamental 
frequency.  * 

Utilizing  the  data  acquired  from  the 
additional  cases  as  described  above,  Fig.  3 
9  plots  the  fluctuating  pressure 
normalized  by  the  dynamic  pressure 
based  on  the  centerline  velocity  UCL  vs. 

L(l  / h ,  the  ratio  of  the  stroke  length  to 
the  orifice  height.  Before  commenting 
on  these  data,  note  that  the  stroke  length 
L0  requires  knowledge  of  the  velocity  -0.4  -0.2  0  0.2  0.4  0.6 

profile.  Since  the  velocity  profiles  are  r/,d 

not  available  for  all  cases,  U-  is  Figure  10.  Normalized  axial  velocity  profiles  at  four  phases  during  a 
estimated  by  computing  the  ratio  cycle  for  S=8,  (o)  Re .  =  45.5 ,  (□)  Re .  =  94 ,  (-)  solution  for  pipe  flow 

U ' j /U CL  f°r  each  Stokes  number  based  with  an  oscillating  A p . 
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on  the  profiles  obtained  for  Cases  1-4.  As  shown  in  Fig.  10,  this  ratio  is  approximately  constant  for  a  given  Stokes 
number  and  invariant  as  the  Reynolds  number  changes.  However,  the  measured  exit  velocity  profile  shape  does  not 
agree  with  the  theoretical  profile  for  oscillatory  pipe  flow30  as  shown  in  Fig.  10,  especially  during  the  suction  stroke. 

Returning  to  Fig.  9,  the  pressure  loss  data  asymptote  to  a  value  of  approximately  0.75  as  Z0//i  increases  beyond 
a  value  of  ~4.  This  suggests  that  when  the  fluid  particle  excursion  or  stroke  length  is  much  larger  than  the  orifice 
height  h ,  minor  “nonlinear”  losses  due  to  entrance  and  exit  effects  dominate  the  flow.  However,  the  magnitude  of 
these  losses  and  the  degree  of  nonlinear  distortion  is  likely  to  be  strongly  dependent  on  Reynolds  number.  For 
shorter  stroke  lengths,  as  in  Case  3  shown  in  Fig.  8(c),  when  the  particle  excursion  is  smaller  than  h ,  the  flow 
presumably  approaches  a  fully-developed  state,  and  major  losses  which  scale  with  u  dominate  the  minor  losses. 
Here,  “fully  developed”  means  that  there  exists  a  region  within  the  orifice  away  from  either  exit,  where  the  velocity 
profile  at  a  given  point  during  the  cycle  is  not  a  function  of  axial  position  x  . 

The  results  of  the  CFD  simulations  allow  us  to  explore  this  hypothesis  further.  Fig.  11  shows  the  variation  of  the 
spanwise  vorticity  for  the  three  computational  cases  (5,  6  ,and  7).  For  Cases  5  &  6,  no  jet  is  formed  as  seen  in  Fig. 
11(a)  and  (b).  The  vorticity  contours  show  that  vortices  formed  during  the  expulsion  cycle  for  these  cases  are 
ingested  back  during  the  suction  cycle,  leading  to  the  trapping  of  vortices  inside  the  orifice  and  increased  losses  over 
and  above  that  predicted  by  linear  theory.  In  contrast,  clear  jet  formation  occurs  for  Case  7  as  shown  in  Fig.  11(c). 
These  results  are  consistent  with  the  jet  formation  criterion  proposed  by  Utturkar28  and  Utturkar  et  al.29,  which  have 
shown  that  jets  are  formed  above  a  critical  value  of  Re(/.S'2  >  0.16  and  1  for  a  shaip  orifice  and  2-D  slot, 

respectively.  The  flow  is  qualitatively  the  same  for  Cases  5  and  6  despite  the  doubling  of  h  ,  suggesting  that  orifice 
height  plays  a  relatively  minor  role  in  determining  jet  formation. 


Figure  11.  Vorticity  contours  during  the  maximum  expulsion  portion  of  the  cycle  for  (a)  Case  5,  (b)  Case  6, 

and  (c)  Case  7. 

Next,  temporal  variations  of  the  pressure  as  a  function  of  location  along  the  slot  and  volume  flow  rate  obtained 
from  the  numerical  simulations  are  shown  in  Fig.  12.  A  comparison  of  the  pressure  traces  shows  a  significant  phase 
shift  in  the  pressure  (see  probe  1  position)  in  the  jet  formation  case  (Case  7)  with  the  largest  stroke  length  compared 
to  Cases  5  and  6.  Furthermore,  all  of  the  pressure  traces  in  Case  7  show  significant  nonlinear  distortion,  a  behavior 
which  is  similar  to  that  seen  for  the  low  Stokes  number  Cases  1  and  2  in  Fig.  8,  both  of  which  correspond  to  jet 
formation.  The  deviation  of  the  pressure  signal  from  a  sinusoidal  shape  is  a  clear  indication  that  the  periodic 
ejection  and  ingestion  of  strong  vortices  has  a  major  effect  on  the  resistance  characteristics  of  the  orifice.  However, 
for  the  highest  Reynolds  number  (Case  4),  the  nonlinear  losses  are  reduced  in  magnitude  such  that  the  pressure 
distortion  in  Fig.  8(d)  is  not  evident. 
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Figure  12.  (a)  Probe  positions  in  numerical  simulations.  Volume  flow  rate/pressure  at  various  locations 

inside  the  orifice  for  (b)  Case  5,  (c)  Case  6  and  (d)  Case  7  during  a  cycle.  The  pressure  at  probes  2  to  6  are 

spatially  averaged  pressure  across  the  orifice. 

The  CFD  data  permits  the  examination  of  the  evolution  of  the  velocity  profile  inside  the  slot.  Fig.  13  shows  the 
computed  velocity  profiles  at  various  locations  in  the  orifice  corresponding  to  the  maximum  expulsion  phase.  As 
with  the  pressure.  Cases  5  and  6  are  qualitatively  similar,  and  both  show  that  the  velocity  profile  undergoes 
significant  development  along  the  length  of  the  slot.  A  close  examination  of  flow  animations  reveals  the  strong 
periodic  influence  of  the  vorticity  as  it  is  expelled  and  then  ingested  again.  This  effect  is  presumably  enhanced  by 
the  sharp  corners  of  the  slot  compared  to  the  rounded  edges  studied  by  Smith  and  Swift.22  In  contrast,  the  evolution 
of  the  velocity  profiles  for  Case  7,  in  which  a  jet  is  formed,  is  significantly  different  from  the  previous  two  cases. 
For  this  case,  where  the  stroke  length  is  much  larger  than  the  orifice  height,  the  core  flow  moves  in  phase  at  each 
location  inside  the  orifice.  Furthermore,  it  appears  that  the  flow  approaches  a  quasi-steady  behavior  for  large  L0/h 
and  is  dominated  by  entrance  and  exit  losses,  the  magnitude  of  which  are  inversely  proportional  to  Reynolds 
number. 
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Figure  13  .  Velocity  profile  at  different  location  inside  the  orifice,  during  the  maximum  expulsion  portion  of 
the  cycle,  (a)  Case  5,  (b)  Case  6  and  (c)  Case  7.  The  locations  are  indicated  in  Figure  12(a). 

VII.  Conclusions  &  Future  Work 

A  thorough  understanding  of  the  unsteady  flow  behavior  and  associated  losses  in  the  orifice/slot  of  zero-net 
mass-flux  devices  is  critical  for  the  development  of  accurate  low-order  models  suitable  for  design.  This 
experimental  and  computational  study  has  examined  ZNMF  devices  exhausting  into  a  quiescent  medium.  The 
geometry  is  characterized  by  a  sharp-edged,  axisymmetric  orifice  or  its  2-D  slot  counterpart.  It  has  been  found  that 
the  flow  field  emanating  from  the  orifice/slot  is  characterized  by  both  linear  and  nonlinear  losses.  Due  to  the  sharp- 
edged  geometry,  the  linear  losses  are  under  predicted  by  the  theoretical  linear  solution  for  pipe/channel  flow  driven 
by  an  oscillating  pressure  gradient.  Furthermore,  the  flow  field  inside  the  orifice  may  vary  significantly  from  the 
linear  solution  due  to  the  persistent  presence  of  vorticity  in  the  orifice  from  previous  expulsion/ingestion  cycles. 
Nevertheless,  it  was  found  that  the  pressure  drop  in  the  orifice,  when  normalized  by  the  dynamic  pressure  based  on 
the  amplitude  of  the  flow  centerline  velocity,  collapsed  at  high  dimensionless  stoke  lengths  over  a  range  of  Reynolds 
and  Stokes  numbers  in  the  experiments  and  simulations.  When  the  stroke  length  is  large  compared  to  the  orifice 
height,  this  dimensionless  pressure  drop  is  small,  0(0.75),  presumably  because  the  flow  is  dominated  by  nonlinear 
minor  losses  associated  with  entrance  effects.  However,  when  the  stroke  length  is  small,  the  pressure  drop  rises, 
presumably  due  to  the  lower  Reynolds  numbers  and  the  major  losses  typical  of  fully-developed  flows. 

In  future  experiments  and  simulations,  additional  Stokes  numbers  and  values  of  h/ d  will  be  investigated.  It  is 
anticipated  that  the  additional  experimental  and  computational  data  will  permit  a  physics-based  empirical  correlation 
for  the  nonlinear  loss  coefficient  as  a  function  of  Reynolds  number  and  stroke  length. 
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